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INTRODUCTION 


The field of sensitivity derivative analysis is emerging as one of the more 
fruitful areas of engineering research. The reason for this is the recognition 
of the many practical uses for sensitivity derivatives. Beyond the historical 
use of derivatives in connection with formal mathematical optimization 
techniques, recent work has been reported in using sensitivity derivatives in 
approximate analysis, assessing design trends, analytical model improvement, and 
determining effects of parameter uncertainties (refs. 1 through 7). 

Work supported by the NASA Langley Research Center, under a grant in sensitivity 
analysis, has been focused on derivatives of thermal response of structures 
(refs. 8 and 9). Most recently, in-house implementations of generalized struc- 
tural sensitivity capability in the SPAR and EAL computer programs (refs. 10 
and 11) have been completed. Work in the sensitivity area is being expanded, 
and recent developments both in and outside the structures area have been sur- 
veyed to guide the future effort. This paper reviews some innovative techniques 
applicable to sensitivity analysis of discretized structural systems. These 
techniques include a finite-difference step-size selection algorithm, a method 
for derivatives of iterative solutions, a Green's function technique for deriva- 
tives of transient response, a simultaneous calculation of temperatures and 
their derivatives, derivatives with respect to shape, and derivatives of optimum 
designs with respect to problem parameters. Computerized implementations of 
sensitivity analysis and applications of sensitivity derivatives are also dis- 
cussed. Finally, some of the critical needs in the structural sensitivity area 
are indicated along with Langley plans for dealing with some of these needs. 


DISCIPLINES CONTRIBUTING TO SENSITIVITY ANALYSIS DEVELOPMENT 


Sensitivity methodology has been and continues to be an important research area 
for many disciplines. Appreciation for the uses of sensitivity analysis by a 
broad spectrum of researchers outside the structures area is very evident. Some 
of those disciplines are indicated in figure 1. For the most part, the 
motivation in these other disciplines is the need to quantify the effect of 
uncertainties in parameters of a system model on the predictions of the model. 
Examples from physical chemistry are described in references 7 and 12 through 18. 
A specific use is given in figure 2. Electronics and control theory represented 
the origin of this type of sensitivity work (refs. 19 and 20) in addition to 
the use of derivatives to synthesize systems. Recent work in physiology with 
both human and bacteriological system models is described in reference 21. In 
the thermodynamics area, reference 22 describes the calculations and use of 
derivatives of the chemical composition with respect to thermodynamic properties 
in the mathematical modeling of a coal gasification process. Finally, analytical 
techniques are beginning to emerge to calculate derivatives of aerodynamic quan- 
tities with respect to flow parameters (refs. 23 and 24) as described by Bristow 
(ref. 25). This paper focuses: on contributions to sensitivity methodology ori- 
ginating in or applicable to the structural analysis field. 


• CHEMICAL KINETICS 

• ELECTRONICS AND CONTROL 

• PHYSIOLOGY 

• THERMODYNAMICS 

• AERODYNAMICS 

• STRUCTURAL ANALYSIS 


Figure 1 


APPLICATION OF SENSITIVITY DERIVATIVES TO ATMOSPHERIC POLLUTION MODEL 


Sensitivity analysis has been used to assess the effects of uncertainties in 
emission and meteorological parameters on the predictions from a mathematical 
model for photochemical air pollution (ref. 7). The atmospheric diffusion 
equation which governs the degree of pollution of a volume of air (an air 
parcel) contains several parameters: (1) mixing depth - the vertical height of 

the air parcel containing pollutants; (2) initial concentration of pollutant; 

(3) photolysis intensity - the rate of photochemical activity; (4) emission 
rate - the rate at which the pollutant is emitted into the air parcel; and (5) 
ambient temperature of the air . The calculation of derivatives of concentrations 
of various pollutants with respect to the aforementioned parameters is described 
in reference 7. The derivatives were used to rank the importance of the parameters. 
The calculations were carried out for the example of an "oxidant episode" which 
occurred in Southern California in 1974. The mathematical simulation of the 
event began in downtown Los Angeles at 5 a.m. and terminated in San Bernardino 
County at 8 p.m. The graph in figure 2 shows the sensitivity of the concentration 
of the pollutant nitric oxide (NO) with respect to each parameter as a function 
of time. Results indicate that early in the episode the initial concentration 
and mixing depth are the most important parameters. Midway through, emission 
rate and ambient temperature were most important, and late in the calculation, 
ambient temperature and photolysis intensity were most critical. These types 
of data indicated the need for more exact measurements of the key parameters to 
improve the air quality mathematical model. 

PARAMETER 



TIME, min 

EXAMPLE PROBLEM - 

POLLUTION DUE TO OXIDANT EPISODE SENSITIVITY OF NO CONCENTRATION 


• STUDY RESULTS LED TO RESEARCH FOR MORE 
PRECISE MEASUREMENTS OF KEY PARAMETERS 

(REF. 7) 


Figure 2 
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OPTIMUM STEP SIZE FOR FINITE-DIFFERENCE DERIVATIVES 


The most straightforward method of calculating derivatives is to use a finite- 
difference approximation. One of the most serious shortcomings of the finite- 
difference method is the uncertainty in the choice of a perturbation step size. 

If the step size is too large, truncation errors may occur. These can be 
thought of as errors due to retention of only the lowest order terms of a Taylor 
series representation of a perturbed function. If the step size is too small, 
condition errors may occur (ref. 26). These errors are due to subtraction of 
nearly equal numbers. In a recent paper (ref. 27), an algorithm was developed 
to determine the optimum finite-difference step-size, i.e., one that balances 
the truncation and condition errors. The algorithm is based on approximating 
the truncation error as a linear function of step size h and the condition 
error as a linear function of 1/h. The optimum step size is obtained by equating 
the condition and truncation errors (fig. 3). This technique has been tested 
on functions which could be differentiated analytically (ref. 27) and was found 
to be very effective. A logical extension of this work would be to apply it to 
matrix equations. 


• WANT BEST ESTIMATE OF ~ ( f ( v + h ) - f(v)) 

IF h TOO LARGE - TRUNCATION ERROR = T( h ) 

IF h TOO SMALL - CONDITION ERROR = C(h) 

EXPRESS T(h ) AND C(h) AS SIMPLE COMPUTABLE FUNCTIONS 
CHOOSE "OPTIMUM" STEP SIZE h SO THAT 
C(h) = T(h ) 

• RESULT - FORMULA FOR h 

• FORMULA VERIFIED BY TESTS ON ANALYTICAL FUNCTIONS 

• NEED TO IMPLEMENT FOR MATRIX EQUATIONS 

(REF. 26) 


THE 

PROBLEM 


THE 

SOLUTION 


Figure 3 


DERIVATIVES OF ITERATIVE SOLUTIONS 


In many structural design problems, the response U is the solution of an 
algebraic system f(U, v) = 0, where v is a design parameter (fig. 4). When 
the system is solved iteratively, the iterative process is terminated when the 
s oluti on error is reduced below a certain tolerance. To obtain the derivative 
of U with respect to a design parameter by finite differences, the parameter is 


. , , ■ . i ■ J ^ uic par ametcr 

solution process is repeated to obtain U b . The derivative 


is 


then approximated by a finite-difference ratio. The error inherent in this pro- 
cess is due to the termination^ of the_iterative solution process before an exact 
solution is obtained. Thus U and U b , obtained by iteration, are only approxi- 
mations to the corresponding exact solutions U and U h , respectively. Because 
of noise in the solution process, the difference between U and Uh can be 
finite, even for very small values of the perturbation h. In fact, the error 
is most severe when small values of h are required to avoid large truncation 
errors in the derivative. A remedy, which is being developed by the second author 
0 j •^ s J paper ’. 1S t0 deflne a modified perturbed solution Uh which satisfies a 
modified equation whose right-hand side is not zero but is the residual of the 
approximate unperturbed equation. By this construction, US approaches 0 as h 
approaches zero. Then U b replaces U h in the derivative formula. Finally, U 
serves as the first approximation in the iteration process for U^ 


• 0 IS SOLUTION TO f(U,v) = 0 

• U h IS SOLUTION TO f (u h , v + 



APPROXIMATE 

0 AND 0 . 
h 

. OBTAINED BY ITERATION 


dU/dv = l/h(u h - u) 


• ERROR MAY BE LARGE DESPITE SMALL h 

• SOLUTION - DEFINE U b SUCH THAT 

f ( U h’ V+ h ) = f( ^’ v) THEN dU/dv = 1 /h (u # h -0) 
RESIDUAL 


• WHY IT INITIAL GUESS FOR U? IS 0 
WORKS h 

U* APPROACHES U FOR SMALL h 


Figure 4 
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GREEN'S FUNCTION METHOD FOR DERIVATIVES OF TRANSIENT RESPONSE 


This method, which is well known in applications to solutions of nonhomogeneous 
differential equations, has been used extensively by physical chemistry 
researchers (refs. 12 through 18) for calculation of derivatives of response 
quantities governed by systems of first-order nonlinear ordinary differential 
equations such as equation (1) shown in figure 5. Numerous applications have 
been performed for chemical kinetics problems related to air pollution studies. 

As indicated in this figure, the derivative of the response vector Y with 
respect to a parameter a satisfies equation (2). The derivative may be repre- 
sented by an integral expression (eq. 3) involving a kernel K which is the 
Green's function. The Green's function is the solution to the initial value 
problem given by equation (4). Comparison of the effort needed to solve equation 
(2) versus (4) indicates that the Green's function technique is advantageous if 
the number of design variables m exceeds the number of equations n in the 
system. One approach to obtaining K is to solve equation (4) directly using 
an implicit numerical integration technique (ref. 12). An alternate solution 
for the Green's function is to use the Magnus method (ref. 17) whereby K is 
expressed as an exponential function of a matrix which is the time integral of 
the Jacobian matrix J. Because the equation of transient heat transfer is a 
special case of equation (1), the Green's function method is directly applicable 
to sensitivity of transient temperatures. It is planned to pursue this line of 
research at Langley as part of our sensitivity development. 


• IMPLEMENTED BY PHYSICAL CHEMISTRY RESEARCHERS 


• GENERAL 
PROBLEM 


^ = f(Y,a.t) 
dY a 

-rr = JY + f 

dt a a 

J = df/dY 


n EQUATIONS 
mn EQUATIONS 
Y q = dY/da 


( 1 ) 

( 2 ) 


• SOLUTION BY 
GREEN'S 
FUNCTION 


t 

Y a = / K(t.T>f a dT 
0 

(t.T ) - JK(t,T) = 0 


n 2 EQUATIONS 


(3) 


(4) 


• APPLICABLE TO TRANSIENT HEAT TRANSFER 


Figure 5 
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CONCURRENT CALCULATIONS OF TRANSIENT TEMPERATURES AND DERIVATIVES 


Derivatives of transient response such as structural displacements and tem- 
peratures have been computed following the calculation of the response itself 
using analytical techniques (refs. 9, 28, 29) and by the Green's function 
technique (ref. 12). Recently an algorithm for concurrent calculation of 
transient temperatures and their finite-difference derivatives has been 
developed (ref. 9). Figure 6 depicts an application to the transient thermal 
response of an insulated cylindrical shell. Temperatures throughout the shell 
are computed using an implicit numerical integration technique. Along with 
temperatures at each time step, finite-difference sensitivity derivatives are 
calculated with respect to design variables representing insulation thicknesses 
at 10 locations on the shell surface. The key to the success of the method 
is that during each time step, when a nonlinear algebraic equation is solved by 
iteration for the current temperature and the perturbed temperature, the same 
time step is used for both solutions. Further, the unperturbed temperature 
serves as the initial guess in the iteration for the perturbed temperature. The 
numbers in the table are solution times in seconds. The results indicate that 
the timesaving from the concurrent calculation is substantial, and nearly a 
factor of 4 advantage is obtained. 


ERROR 

TOLERANCE 

ON 

TEMPERATURE 

SEQUENTIAL 

CONCURRENT 

CPU TIME 
FOR TEMPERATURE 
CALCULATION 

CPU TIME 
FOR EACH 
DERIVATIVE 

CPU TIME 
FOR TEMPERATURE 
PLUS TEN DERIVATIVES 

CPU TIME 
FOR EACH 
DERIVATIVE 

.001 

113 

113 

500 

39 

.003 

78 

78 

296 

22 

.010 

46 

46 

168 

12 



TYPICAL DERIVATIVE 
dT 


T = TEMPERATURE 

t. = INSULATION THICKNESS 
ins 


Figure 6 
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SENSITIVITY DERIVATIVES FOR SHAPE DESIGN VARIABLES 


A relatively new topic in structural sensitivity analysis is the calculation of 
derivatives with respect to shape design variables. Examples are derivatives of 
displacements or stresses with respect to a beam length or a membrane area 
(fig. 7). Two approaches have been used. One approach is to differentiate the 
discretized equations resulting from a finite-element representation. A draw- 
back to this technique is that when shape design variables change, the finite- 
element mesh is modified. The resulting mesh distortion changes the discretization 
* error and leads to inaccurate derivatives. The second approach, which avoids mesh 

► distortion errors, is to reverse the order of differentiation and discretization 

(refs. 30 through 32). The procedure is to differentiate the continuum equations 
of equilibrium and discretize the resulting integral equations. This method uses 
the concept of a material derivative from continuum mechanics which is composed 
of two parts: a derivative corresponding to a fixed shape, and a contribution 

from the change of the boundary. The preferred choice between the two methods 
is not yet clear. The second approach avoids mesh distortion by its formulation 
but does not permit shape differentiation of a discretized set of equations. 

The first method, although suffering from. mesh distortion errors, could benefit 
from a built-in adaptive mesh generation capability which would reduce the mesh 
distortion. 


• EXAMPLES - DERIVATIVE OF 

•DISPLACEMENTS WITH LENGTH OF BEAM 

RESPECT 

• STRESSES T0 AREA OF MEMBRANE 


• FIRST METHOD - DISCRETIZE FIRST THEN DIFFERENTIATE 

• NUMERICAL ERRORS DUE TO MESH DISTORTION 

• REDUCE ERRORS BY ADAPTIVE MESH GENERATION 


• SECOND METHOD - BASED ON MATERIAL DERIVATIVE 

• DIFFERENTIATE CONTINUUM EQUATIONS THEN DISCRETIZE 
•AVOIDS MESH DISTORTION ERRORS 


Figure 7 
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SENSITIVITY OF OPTIMUM DESIGNS TO PROBLEM PARAMETERS 


The problem addressed by this technique (refs. 33, 34) is to obtain derivatives 
of an objective function F and design variables V from an optimized solution 
with respect to parameters P which were held constant during the optimization. 
The most obvious and thus most useful application of the technique is extra- 
polation of an optimum design for variations of a problem parameter. For example, 
the effect of varying the height H of the truss in figure 8 is assessed by using 
optimum sensitivity derivatives. Extrapolated values of the mass F and one 
of the cross-sectional areas A., based on derivative with respect to H, are 
compared with those obtained byreoptimization with different values of H. As 
shown in the lower right portion of figure 8, the results agree very closely for 
up to a 20-percent change in H. Other applications of these types of derivatives 
include optimization for multiple-objective functions, assessing the effects of 
adding or deleting constraints, and most recently using the derivatives as links 
between subsystems during multilevel optimization (ref. 35). 


GIVEN AN OPTIMUM DESIGN: 

F = OBJECTIVE FUNCTION V = DESIGN VARIABLES 

P = PROBLEM PARAMETERS 


ANALYSIS GIVES -|jj AND ^ 



• EXAMPLE - EFFECT OF TRUSS HEIGHT 
ON OPTIMUM DESIGN 


A- EXTRAPOLATION 



(REF. 33) 


AH/H, percent 


Figure 8 


376 


COMPUTER IMPLEMENTATION 


Some progress has occurred in providing general-purpose software for sensitivity 
analysis (fig. 9). The Green's function technique described earlier has been imple- 
mented in a computer program denoted AIM (ref. 18). Use of this program requires 
supplying subroutines to define the system of equations - specifically, the 
vector f and the matrix J in figure 5. The capability for computing deriva- 
tives of static displacements, stresses, and vibration and buckling eigenvalues 
and eigenvectors has been implemented in the SPAR finite-element program 
(ref. 10) and EAL (refs. 11, 36). The EAL (Engineering Analysis Language) sys- 
tem contains the SPAR finite-element modules, but additionally EAL provides 
FORTRAN-like commands which permit branching, testing data, looping, and calling 
the SPAR modules (similar to calling FORTRAN subroutines). A recent level of 
a proprietary version of the NASTRAN computer program also has capability for 
static displacement, stress, and eigenvalue derivatives (ref. 37). 


• AIM (GREEN'S FUNCTION TECHNIQUE) 

•GENERAL FIRST-ORDER EQUATIONS 

• DERIVATIVES WRT PARAMETERS IN EQUATIONS 

• SPAR (COSMIC) AND EAL (PROPRIETARY) 

•DERIVATIVES OF - 

•DISPLACEMENTS 

•STRESSES 

•EIGENVALUES 

•EIGENVECTORS 

• NASTRAN (PROPRIETARY) 

•DERIVATIVES OF - 
•DISPLACEMENTS 
•STRESSES 
•EIGENVALUES 


Figure 9 
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APPLICATION OF STRUCTURAL SENSITIVITY ANALYSIS TO SPACE ANTENNA 


An application of sensitivity analysis to reveal structural design trends is 
illustrated in figure 10. The structure is an Earth-orbiting antenna 
reflector subjected to nonuniform heating leading to thermal distortions which 
can degrade antenna performance (ref. 11). The structure is modeled using only 
rod elements. There are three design variables representing, respectively, the 
cross-sectional areas of the elements in the upper surface (A 3 ), the elements 
joining the upper and lower surfaces (A 2 ), and the elements in the lower 
surface (A 3 ). Derivatives of the center deflection with respect to each 
design variable were calculated and are shown in the figure. A positive 
derivative indicates that increasing the design variable increases the 
response. A negative derivative indicates that increasing the design variable 
decreases the response. The seemingly contradictory result that increasing a 
design variable can increase a response stems from the fact that the thermal 
loads are proportional to the rod cross-sectional areas. From the table in 
figure 10 , we see that increasing Aj has the largest effect on reducing deflection 
but at the cost of a weight increase. On the other hand, decreasing either An or 
A 3 would reduce the deflection and at the same time reduce weight. It is at the 
discretion of the designer as to which of the alternatives is a better choice. 

The sensitivity derivatives provide the data for that judgment. 



ANTENNA MODEL 

w = CENTER DEFLECTION 
0 


DESIGN VARIABLE - 
ELEMENT AREA, Aj 

dw 0 /SAj 

UPPER SURFACE 

-2.4 x 10" 4 

DIAGONALS 

8.3x 10“ 5 

LOWER SURFACE 

1.8x lO" 4 


DERIVATIVES OF CENTER DEFLECTION 


Figure 10 
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APPLICATION OF SENSITIVITY ANALYSIS TO SHUTTLE PAYLOAD 




4 * 


^orl t1v1ty u? na l ysis has been used t0 rec *esign the Airborne Support Equipment 
(ASE) assembly on the Space Shuttle orbiter (ref. 5). The ASE supports the 
inertial upper stage (IUS) vehicle in the payload bay. The purpose of the IUS 
vehicle is to transport payloads further into space once the orbiter has 
established low-Earth orbit. The shell of the ASE is subjected to large launch 
loads and is designed for an ultimate load of 3000 Ib/in. An initial sculptured 
skin design met the design load but was too heavy. A sensitivity analysis was 
performed to determine which skin gages had the largest effect on loads and to 
determine which type of modified construction would give the largest weight 
reduction (among isogrid, waffle, and stiffened skin). The structure was 
modeled and divided into zones as shown on the right side of figure 11. 
Derivatives of compressive loads and weight with respect to longitudinal, 
circumferential, and shear stiffening were computed for each zone. As an 
example, consider derivatives with respect to the longitudinal stiffness design 
variable t a . The numbers in the zones are ratios of derivatives of load to 
derivatives of weight for the sculptured skin design. Negative values indicate 
that increasing a design variable decreases the load, and positive values 
indicate that increasing a design variable increases the load. The analysis 
revealed that derivatives with respect to longitudinal and shear stiffness were 
the largest, derivatives with respect to circumferential stiffness were 
negligible, and derivatives with respect to shear stiffness were nearly all 
positive. Based on these results, the shell was redesigned as a 
longitudinally stiffened machined skin. (The isogrid was rejected because of 
high shear stiffness; the waffle construction was rejected due to unneeded high 
circumferential stiffness.) The resulting design satisfies the ultimate load 
constraint with a large margin of safety and an acceptably low weight. 



BASELINE ASE CONFIGURATION 



DISTRIBUTION OF DERIVATIVE RATIOS 


• STUDY RESULTS 

•LARGE DERIVATIVES WRT 
LONGITUDINAL AND SHEAR 
STIFFNESS 

• SMALL DERIVATIVES WRT 
CIRCUMFERENTIAL STIFFNESS 


dN/dt a 

a w/dt a 

t a = LONGITUDINAL STIFFNESS 


• OUTCOME - REDESIGNED SHELL AS LONGITUDINALLY STIFFENED MACHINED PANEL 
SATISFIED DESIGN REQUIREMENTS WITH LOW WEIGHT 

(REF. 5) 

Figure 11 
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STRUCTURAL SENSITIVITY ANALYSIS NEEDS 


As a result of surveying methods applicable to computing structural sensitivity 
derivatives, a list of needs has emerged (fig. 12). First, continued development 
of methods for derivatives of transient response and derivatives with respect to 
shape design variables and material properties should have high priority. 

Further, techniques developed for sensitivity derivatives in nonstructural 
disciplines such as physical chemistry have much to offer and should be 
evaluated for their adaptability to structural areas. It appears that 
structural designers have made insufficient use of the power and utility of 
sensitivity derivatives to guide design modifications and to assess 
uncertainties in their models. Their use can be accelerated by demonstrations 
of practical applications of sensitivity analysis and careful documentation (by 
optimization and sensitivity specialists) to guide structural analysts and 
designers not experienced in formal optimization and sensitivity analysis. 
Finally, sensitivity analysis needs to be routinely included as a standard 
feature in general-purpose structural analysis software packages. Near-term 
plans at Langley include evaluation of the Green's function method for 
derivatives of transient thermal response, methods for derivatives of spacecraft 
thermal response with respect to material properties, and implementation of the 
optimum finite-difference step-size technique for finite-element sensitivity 
analysis. Concurrent with this effort, demonstration problems will be selected 
and solved. 

• TRANSIENT RESPONSE 

• MATERIAL PROPERTIES AND SHAPE 

• PRACTICAL APPLICATIONS 

• ROUTINE INCLUSION IN GENERAL-PURPOSE COMPUTER PROGRAMS 


LANGLEY PLANS 

• EVALUATE GREEN'S FUNCTION METHOD FOR TRANSIENT TEMPERATURES 

• DERIVATIVES WRT MATERIAL PROPERTIES 

• IMPLEMENT FINITE-DIFFERENCE STEP-SIZE ALGORITHM 

• SENSITIVITY DEMONSTRATION PROBLEMS 


Figure 12 
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SUMMARY 


This paper was based on a recently conducted survey of methods for sensitivity 
analysis of structural response (fig. 13). The survey was not limited to 
research in the structural area alone and revealed that a broad range of dis- 
ciplines are using sensitivity analysis and contributing to the methodology. 

In almost every instance, methods from the nonstructural disciplines are dir- 
ectly applicable to the structures area. An example application from chemical 
kinetics was described in which sensitivity analysis was used to assess the 
impact of parameter uncertainties in a mathematical model used in air pollution 
studies. 

The bulk of the paper has focused on a selected set of innovative methods 
applicable to sensitivity analysis of structural systems. The analysis 
techniques include a finite-difference step-size selection algorithm, a method 
for derivatives of iterative solutions, a Green's function technique for 
derivatives of transient response, concurrent calculation of temperatures and 
their derivatives, derivatives with respect to shape, and derivatives of optimum 
designs with respect to problem parameters. Two applications were described 
wherein derivatives were used to guide structural design changes to improve an 
engineering design without recourse to formal mathematical optimization. Plans 
at Langley for contributing to identified critical needs were cited. Among the 
needs were implementation of methods for derivatives of transient response, 
derivatives with respect to shape and material properties, solution and 
documentation of sensitivity analysis demonstration problems, and routine 
inclusion of sensitivity analysis as a feature in general-purpose structural 
analysis computer programs. Langley near-term plans in the sensitivity area 
include evaluating the Green's function method for derivatives of transient 
thermal response, developing methods for derivatives with respect to material 
properties, implementation of the finite-difference step-size algorithm, and 
solution of sensitivity demonstration problems. 


• LANGLEY CONDUCTING SURVEY OF METHODS FOR SENSITIVITY DERIVATIVES 

• PAPER REVIEWED RECENTLY DEVELOPED TECHNIQUES 

•OPTIMUM FINITE-DIFFERENCE STEP-SIZE ALGORITHM 
•METHOD FOR DERIVATIVES OF ITERATIVE SOLUTIONS 
•GREEN'S FUNCTION METHOD FOR TRANSIENT RESPONSE 

• SIMULTANEOUS CALCULATION OF TEMPERATURES AND DERIVATIVES 
•DERIVATIVES WITH RESPECT TO SHAPE 

• DERIVATIVES OF OPTIMUM DESIGNS WITH RESPECT TO PROBLEM 
PARAMETERS 

• REVIEWED USES OF DERIVATIVES AND CITED APPLICATIONS 

• CITED NEEDS AND OUTLINED LANGLEY PLANS 

Figure 13 
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